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Abstract 

A re-examination of the energy cosmic ray spectrum above 10 20 eV is presented. 
The overall data-base provides evidence, albeit still statistically limited, that non- 
nucleon primaries could be present at the end of the spectrum. In particular, the 
possible appearance of superheavy nuclei (seldom discussed in the literature) is 
analysed in detail. 
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The origin and nature of cosmic radiation have been a constant source of 
mystery and discovery since 1949 [1]. Most notably, Greisen, Zatsepin and 
Kuz'min (GZK) pointed out that extremely high energy cosmic rays (usually 
assumed to be nucleons or nuclei) undergo reactions with the pervasive mi- 
crowave background radiation (MBR) yielding a steep drop in their energy 
attenuation length [2]. Specifically, any proton energy above 50 EeV is de- 
graded by resonant scattering via 7 + p — > A — > p/n + tt, and heavy nuclei 
with energies above a few tens EeV get attenuated mainly by photodisintegra- 
tion off the MBR and intergalactic infrared background photons (IR). Over 
the last few years, several giant air showers have been detected which confirm 
the arrival of particles with energies > 100 EeV, this is, above the GZK cut- 
off (see [3] for a recent survey). Many models have been proposed as source 
candidates of such high energy events [4] , however, it is not known for certain 
at the present time from where the rays originate. 

In revealing their origin, the observed anisotropy of these cosmic rays is one 
of the most useful features. Very recently, the Fly's Eye [5] and Akeno Giant 
Air Shower Array (AGASA) [6] experiments reported a small but statistically 
significant anisotropy 0(4%) in the cosmic ray flux towards the galactic plane 
at energies around 1 EeV. With increasing energy the picture looks rather 
different: although at E > 40 EeV an enhancement of the flux from the 
Supergalactic plane was reported [7], the arrival directions above 100 EeV are 
best described as isotropic, without imprint of correlation with the galactic 
plane or Supergalactic plane [8]. There are two extreme explanations for this 
puzzle: i) the bunch of nearby sources follows an isotropic distribution (which 
hardly could be the case) ii) One (A few) source(s) dominates at the highest 
energies whilst the background fields of the intergalactic medium strongly 
modify the particle propagation. For the latter explanation, it was suggested 
that a Galactic wind akin the solar wind could bend all the orbits of the 
highest energy cosmic rays towards the Virgo cluster ( VC) [9] . Actually, if one 
assumes that these particles are protons, except for the two highest energy 
events (the one recorded at AGASA [10] and the super-GZK event reported 
by the Fly's Eye group [11]) all trajectories can be traced to within less than 
about 20 degrees from Virgo.Q 

At the highest energies, observed extensive air showers seem to be consistent 
with nucleon primaries but due to the poor statistics and large fluctuations 
from shower to shower an accurate determination of the particle species is not 
possible at the moment. Furthermore, extensive air shower simulations depend 



Notice that the highest energy Yakutsk event was excluded from this sample 
because of the great uncertainty on its energy determination. While first estimates 
suggested a primary energy around 120 EeV [12], a re-estimation of the number of 
charged particles at 600 m from the shower core yields a possible primary energy of 
300 EeV [13]. 
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to some extent on the hadronic interaction event generator which complicates 
the interpretation of data even more [14]. Interestingly enough, however, the 
muon component of the highest AGASA event agrees with the expectation 
extrapolated from lower energies [10]. Indeed, a population of piled-up protons 
is expected at 50 EeV [15], and the picture seems quite consistent. On the 
other hand, the Fly's Eye event occurs high in the atmosphere, and, although 
a primary proton cannot be excluded, a heavy nucleus more closely fits its 
shower development [16]. 

It is widely believed that the cosmic ray spectrum beyond the "crossover en- 
ergy" (energy at which the local spectrum becomes comparable to or less 
than the cosmological component) could be associated with the presence of a 
particularly bright extragalactic, though relative nearby source, superimposed 
on a cosmological diffuse background. In Fig. 1 we show the evolved energy 
spectrum of nucleons assuming a cosmologically homogeneus population of 
sources - usually referred to as the universal hypothesis (UH)- [17], together 
with a compilation of recent air shower data [18]. In addition, we show the 
modified spectrum for the case of an extended source described by a Gaussian 
distribution of width 2 Mpc at a distance of 18.3 Mpc (see [19] for details). 
Assuming that there is no other significant energy loss mechanism beyond in- 
teractions with the MBR for cosmic rays traversing parts of the cluster, this 
could be taken as a very crude model of Virgo. It is important to stress that 
for a galactic magnetic field B ga x = 7/xG (as in [9]) which extends -Rhaio ~ 1-5 
Mpc in the galactic halo, the mean flight time of the protons during their trip 
through the Milky Way is ~ 5.05 x 10 6 yr [21]. This means that the bending 
does not add substantially to the travel time, and the continuous energy loss 
within the straight line approximation is expected to be reasonable for the 
problem at hand. From Fig. 1 we realize that the spectrum of the VC suc- 
cessfully reproduces the AGASA data above 100 EeV. However, it apparently 
cannot account for the super-GZK Fly's Eye event. The interpretation that 
we give for this result is that, without specific knowledge of the chemical com- 
position, the best guess is that at the end of the spectrum two different types 
of characters are playing.0 

At this stage, it is interesting to note that the measured density profile of 
the highest energy Yakutsk event (excluded in the analysis of [9]) shows a 
huge number of muons. Remarkably, its arrival direction coincides with the 
300 EeV Fly's Eye event, within angular resolution, possibly indicating a com- 
mon origin. If this is the case, the almost completely muonic nature of this 
event, recently associated with a dust grain impact [23], could be, perhaps, 
the signature of a super-heavy nucleus. 



2 We remark that AGASA data could be also reproduced if sources of ultra high 
energy protons trace the inhomogeneous distribution of luminous matter in the local 
present-epoch universe [22]. 
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It has been generally thought that 56 Fe is a significant end product of stel- 
lar evolution and higher mass nuclei are rare in the cosmic radiation. Strictly 
speaking, the atomic abundances of middle-weight (60 < A < 100) and heavy- 
weight (A > 100) elements are approximately 3 and 5 orders of magnitude 
lower, respectively, than that of the iron group [24]. The synthesis of the 
stable super-heavy nuclides is classically ascribed to three different stellar 
mechanisms referred to as the s-, r-, and p-processes. The s-process results 
from the production of neutrons and their capture by pre-existing seed nu- 
clei on time scales longer than most /3-decay lifetimes. There is observational 
evidence that such a kind of process is presently at work in a variety of chem- 
ically peculiar Red Giants [25] and in special objects like FG Sagittae [26] 
or SN1987A [27]. The abundance of well developed nuclides peaks at mass 
numbers A = 138 and A = 208. The neutron-rich (or r- nuclides) are synthe- 
sized when seed nuclei are subjected to a very intense neutron flux so that 
/9-decays near the line of stability are far too slow to compete with the neu- 
tron capture. It has long been thought that appropriate r-process conditions 
could be found in the hot (T > 10 10 K) and dense (p ~ 10 10 — 10 11 g/cm 3 ) 
neutron-rich (neutronized) material located behind the outgoing shock in a 
type II supernova event [28]. Its abundance distribution peaks at A = 130 
and A = 195. The neutron-deficient (or p-nuclides) are 100 - 1000 times less 
abundant than the corresponding more neutron rich isobars, while their dis- 
tribution roughly parallels the s- and r- nuclides abundance curve. It is quite 
clear that these nuclides cannot be made by neutron capture processes. It is 
generally believed that they are produced from existing seed nuclei of the s- 
or r-type by addition of protons (radiative proton captures) , or by removal of 
neutrons (neutron photodisintegration). The explosion of the H-rich envelopes 
of type II supernovae has long been held responsible for the synthesis of these 
nuclides [24]. 

In light of the above, starbursts appear (hopefully) as the natural sources able 
to produce relativistic super-heavy nuclei. These astrophysical environments 
are supposed to comprise a considerable population of O and Red Giant stars 
[29], and we believe the supernovae rate is as high as 0.2 - 0.3 yr _1 [30]. Of 
special interest here, the arrival directions of the Fly's Eye and Yakutsk super- 
GZK events {b = 9.6°, I = 163° and 6 = 3°,/ = 162°) seem to point towards 
the nearby metally-rich galaxy M82 (b = 41°, I = 141°) [31] which has been 
described as the archetypal starburst galaxy [32] and as a prototype of super- 
wind galaxies [33]. The joint appearance of the galactic wind and the galactic 
magnetic field during particle propagation could certainly account for the re- 
quired 37° deflection. In addition, it was recently suggested that within this 
type of galaxies, iron nuclei can be accelerated to extremely high energies if a 
two step process is invoked [34]. In a first stage, ions are diffusively accelerated 
up to a few PeV at single supernova shock waves in the nuclear region of the 
galaxy [35]. Since the cosmic ray outflow is convection dominated, the typical 
residence time of the nuclei in the starburst results in t ~ 1 x 10 11 s. Thus, 
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the total path traveled is substantially shorter than the mean free path (which 
scales as A~ 2 ^ 3 ) of a super- heavy nucleus (for details see [34]). Those which are 
able to escape from the central region without suffering catastrophic interac- 
tions could be eventually re-accelerated to superhigh energies at the terminal 
shocks of galactic superwinds generated by the starburst. The mechanism ef- 
ficiently improves as the charge number Z of the particle is increased. For 
this second step in the acceleration process, the photon field energy density 
drops to values of the order of the cosmic background radiation (we are now 
far from the starburst region). The dominant mechanism for energy losses in 
the bath of the universal cosmic radiation is the photodisintegration process 
[36]. Notice that the energy loss rate due to photopair production could be 
estimated as Z 2 /A times higher than that of a proton with the same Lorentz 
factor, and thus could be safely neglected [37]. The disintegration rate R (in 
the system of reference where the MBR is at 2.73K) of an extremely high 
energy nucleus with Lorentz factor T, propagating through an isotropic soft 
photon background reads [38] ,Q 





where a stands for the total photon absortion cross section. The density of 
the soft photon background n(e) can be modeled as the sum of: i) the MBR 
component which follows a Planckian distribution of temperature ~ 2.73K, 
ii) the IR background photons as estimated in [39] , iii) a black body spectrum 
with T = 50007^ and a dilution factor of 1.2 x 10~ 15 to account for the optical 
(O) photons. The total photon absortion cross section is characterized by a 
broad maximum, designated as the giant resonance, located at an energy of 
12-20 MeV depending on the nucleus under consideration. For the medium and 
heavy nuclei, A > 50, the cross section can be well represented by a single, 
or in the case of the deformed nuclei, by the superposition of two Lorentzian 
curves of the form 

e' 2 r 2 

a(e ' } = °\ e l-^ ? l e a Tl (2) 



In order to make some estimates, hereafter we refer our calculations to a gold 
nucleus (the resonance parameters are listed in table I [40]). In Fig. 2 we 
show the 197 Au photodisintegration rate due to interactions with the starlight 
and relic photons. At the highest energies, the energy losses are dominated 
by collisions with the tail of 2.73-fT Planckian spectrum. It is straightforward 
to show that a superheavy nucleus of a few hundred EeV emitted by M82 

3 Primed quantities refer to the rest frame of the nucleus. 
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Table 1 



Giant dipole resonance parameters 



e [MeV] cr [mb] 


T [MeV] 


13.15 255 


2.9 


13.90 365 


4.0 



can traverse almost unscathed through the primeval radiation to produce an 
extensive air shower after interaction with the earth atmosphere. 

Additional support for the superheavy nucleus hypothesis comes from the 
CASA-MIA experiment [41] (See in particular Fig. 9). The collected cosmic 
ray data between 10 14 - 10 16 eV tends to favor a supernova shock wave accel- 
eration scenario. The average mass increases with energy, becoming heavier 
above 10 15 eV. At the maximum energy the results are consistent at la level 
with nuclei heavier than iron. However, "lore" has settled down some compar- 
isons of the admittedly limited ultra high energy cosmic ray sample against 
hadronic-interaction-event generators which predicts the arrival of particle 
species heavier than iron [41,42]. We would like to stress that since simula- 
tions are used to interpret data, and then the data is used to modify the 
simulation, one has to be very careful and as we have shown, it is by no means 
clear that superheavy nuclei could not be present at the end of the spectrum. 

The energy spectrum of nearby (around 3Mpc) nuclear sources was discussed 
elsewhere [43] . The analysis showed that particles tend to pile up between 240 
- 270 EeV. This bump-like feature is followed by a simultaneous drop in the 
cosmic ray flux of the preceding bins of energy, changing the relative detection 
probabilities. As a consequence particles in the pile-up are 50% more probable 
than those at lower energies. 

In summary, the recently reported AGASA data can be successfully repro- 
duced by a power law spectrum of nucleons hailing from the VC superim- 
posed on a cosmological diffuse background. The Fly's Eye observations may 
also fit this scenario, albeit with large errors. One might also consider less 
likely astrophysical sources. In particular, our analysis seems to indicate that 
the next-door galaxy M82 could be responsible for some events at the end of 
the CR spectrum. This has also been suggested elsewhere [31,34,43]. At least 
some of these super-GZK events could be due to heavy, and even superheavy 
nuclei. Clearly, more data is needed before this hypothesis can be verified. In 
this regard, the coming avalanche of high quality cosmic ray observations at 
the Southern Auger Observatory [44] will provide new insights to the ideas 
discussed in this letter. 

Note added: After we finished this work, it was argued that the Galactic wind 
model assumed in Ref. [9] is alone responsible for the focusing of positive 
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particles towards the North galactic pole. Therefore the apparent clustering 
of the back-traced CR cannot be interpreted as evidence for a point source, 
this point source identified as M87 [45]. It should be pointed out that the main 
input parameters for the determination of the CR-spectrum in Fig. 1 are the 
spectral index of the source, and the propagation distance of the nucleons in 
the extragalacic medium. The Galactic wind model is just used to collect all 
the traces in only one single direction in the sky. Therefore the discussion 
presented in this letter strongly supports older suspicions regarding M87, like 
the model proposed in Ref. [9]. 
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Fig. 1. The cosmic ray flux spectrum derived from AGASA (square) and Fly's 
Eye (triangle) experiments shown with the shape of the universal hypothesis (UH) 
spectrum (spectral index 7 = 3.27 [20]). We also show the expected flux of ultra 
high energy nucleons from the Virgo cluster (VC). 



197 a 

Au 





















10 
E [eVJ 



Fig. 2. Fractional energy loss for Au photodisintegration on MBR, IR, O, as well 
as the total (solid line). 
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